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Abstract 

Many insects rely on bacterial symbionts with tiny genomes specialized for provisioning nutrients lacking in host diets. Xylem sap and 
phloem sap are both deficient as insect diets, but differ dramatically in nutrient content, potentially affecting symbiont genome 
evolution. For sap-feeding insects, sequenced symbiont genomes are available only for phloem-feeding examples from the suborder 
Sternorrhyncha and xylem-feeding examples from the suborder Auchenorrhyncha, confounding comparisons. We sequenced ge- 
nomes of the obligate symbionts, Sulcia muelleri and Nasuia deltocephalinicola, of the phloem-feeding pest insect, Macrosteles 
quadrilineatus (Auchenorrhyncha: Cicadellidae). Our results reveal that A/asu/a-ALF has the smallest bacterial genome yet sequenced 
(1 1 2 kb), and that the Sulcia-ALF genome (1 90 kb) is smaller than that of Sulcia in other insect lineages. Together, these symbionts 
retain the capability to synthesize the 10 essential amino acids, as observed for several symbiont pairs from xylem-feeding 
Auchenorrhyncha. Nasuia retains genes enabling synthesis of two amino acids, DNA replication, transcription, and translation. 
Both symbionts have lost genes underlying ATP synthesis through oxidative phosphorylation, possibly as a consequence of the 
enriched sugar content of phloem. Shared genomic features, including reassignment of the UGA codon from Stop to tryptophan, 
and phylogenetic results suggest that Nasuia- ALF is most closely related to Zinderia, the betaproteobacterial symbiont of spittlebugs. 
Thus, Nasuia/Zinderia and Sulcia likely represent ancient associates that have co-resided in hosts since the divergence of leafhoppers 
and spittlebugs >200 Ma, and possibly since the origin of the Auchenorrhyncha, >260 Ma. 
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Introduction 

Obligate symbioses with bacteria have enabled insects to ex- 
ploit novel niches, previously out of adaptive reach (Moran 
2007). Prominent examples include plant sap-feeding insects 
in the order Hemiptera that have evolved partnerships with a 
diverse consortium of bacterial symbionts (Buchner 1965). By 
providing essential nutrients lacking in xylem sap or phloem 
sap, these relationships have permitted insect hosts to tap into 
these specialized diets (Buchner 1 965; Baumann 2005; Moran 
et al. 2008). Symbiont acquisition has thereby enabled the 
Hemiptera to diversify into a species-rich and geographically 
ubiquitous insect order. 

Specialization on xylem or phloem in sap-feeding Hemi- 
ptera is a conserved trait, corresponding to ancient clades 
(Grimaldi and Engel 2005; Cryan and Urban 2012). Both 
diet types are deficient in most essential amino acids and 



many cofactors, but otherwise they have vastly different nu- 
tritive qualities (Mattson 1980; Andersen et al. 1989; Sand- 
strom and Petterson 1994; Sandstrom and Moran 1999). In 
general, phloem is more nutritious, with 10-fold greater levels 
of available organic nitrogen and a 1000-fold greater level of 
carbohydrates, mostly sugars (Mattson 1980; Andersen et al. 
1989). However, the influence of phloem versus xylem sap 
diets on the evolution of symbioses remains unknown. 
Because phloem-feeding systems have been studied only in 
the hemipteran suborder, Sternorrhyncha (e.g., Buchnera- 
aphid symbiosis; Shigenobu and Wilson 2011), and xylem- 
feeding examples only in the Cicadomorpha (McCutcheon 
and Moran 2010), which is one of the two clades comprising 
the hemipteran suborder Auchenorrhyncha, diet type and 
phylogenetic relationships are confounded in comparisons. 

Species in the Auchenorrhyncha rely on two obligate sym- 
bionts: Sulcia muelleri (Bacteroidetes), which has been present 
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within Auchenorrhyncha since it emerged 260-280 Ma 
(Moran et al. 2005), and an array of co-primary symbionts 
from different bacterial divisions (fig. 1). In several groups of 
xylem feeders, genomes of Sulcia and its partner symbiont 
have been sequenced; these include cicadas (Sulcia plus Hodg- 
kinia cicadicola [Alphaproteobacteria]; McCutcheon et al. 
2009b), spittlebugs {Sulcia plus Zinderia insecticola [Betapro- 
teobacteria]; McCutcheon and Moran 2010; Koga et al. 
2013), and sharpshooter leafhoppers (Sulcia plus Baumannia 
cicadellinicola [Gammaproteobacteria]; Wu et al. 2006). In 
each case, Sulcia can produce most of the 10 essential 
amino acids, but is unable to produce histidine or methionine. 
And, in each case, these missing pathways are perfectly com- 
plemented by those encoded in the co-primary symbionts (re- 
viewed by McCutcheon and Moran 2010). This is in contrast 
to the aphid-£i/c/7nera model system in which the provisioning 
of all 10 essential amino acids is accomplished by a single 
obligate symbiont (Shigenobu et al. 2000, reviewed by Bau- 
mann 2005). Like aphids, most groups of Auchenorrhyncha 
are phloem feeders, including the treehoppers and most leaf- 
hopper groups (Membracoidea: Membracidae, most Cicadel- 
lidae), and most planthopper groups (e.g., Fulgoroidea: 
Cixiidae, Delphacidae, Flatidae, Fulgoridae) (Dietrich 2003). 
Most of these are also host to dual bacterial symbionts (fig. 
1), based on microscopy and 16S rRNA screening (Buchner 
1965; Noda et al. 2012; Urban and Cryan 2012; Koga et al. 
2013), but to date none has been studied using genomic 
approaches. 

A central aspect of symbiont evolution is extreme genomic 
reduction resulting from maternal transmission and the intern- 
ment of symbionts within specialized host cells (bacteriocytes) 
for millions of years. Repeatedly, symbionts lose pathways and 
genes involved in every functional category, and particularly in 
functions related to cellular regulation, DNA repair, cell enve- 
lope biosynthesis, and secondary metabolism (Shigenobu 
et al. 2000; Tamames et al. 2007; Moran et al. 2008; 
McCutcheon and Moran 2010, 2012). Indeed, auchenor- 
rhynchan symbionts possess the smallest known cellular ge- 
nomes. In some cases, these genomes appear insufficient to 
support cellular life (e.g., McCutcheon and von Dohlen 201 1 ; 
McCutcheon and Moran 2012), yet they retain a subset of 
genes for synthesizing essential amino acids required by hosts 
(McCutcheon and Moran 2012). 

The similarity in amino acid biosynthetic abilities of different 
Sulcia genomes suggests that Sulcia lost pathways for biosyn- 
thesis of methionine and histidine early during the diversifica- 
tion of the Auchenorrhyncha. Sulcia has continuously 
maintained associations with symbionts providing these func- 
tions, even as these partners were sometimes replaced over 
evolutionary time in particular lineages (fig. 1 ; Buchner 1965; 
Moran et al. 2005; McCutcheon and Moran 201 0). It remains 
unclear whether the original partner of Sulcia persists in any 
modern Auchenorrhyncha. However, Zinderia of spittlebugs, 
Nasuia deltocephalinicola of deltocephaline leafhoppers (Noda 



et al. 2012), and Vidania fulgoroideae of fulgorid planthop- 
pers (Gonella et al. 201 1) form a clade within Betaproteobac- 
teria in 16S rDNA-based phylogenies. This BetaSymb lineage 
has been proposed to represent the ancestral partner of Sulcia 
that has been replaced in several descendant lineages of the 
Auchenorrhyncha, perhaps in response to dietary or ecological 
shifts by hosts (see fig. 1 ; Koga et al. 201 3). Conclusions from 
sequence-based phylogenies are hampered by extreme accel- 
eration in sequence evolution and extreme nucleotide base 
compositional biases of symbionts. Shared genomic features 
across Zinderia, Nasuia, and Vidania could provide more de- 
finitive support for BetaSymb as a clade and as an ancestral 
symbiont of Auchenorrhyncha. But, to date, only a single 
Zinderia genome has been available. 

To address the questions of 1) the pattern of genome evo- 
lution in dual symbionts of a phloem-feeding host, and 2) the 
origin for the BetaSymb lineage, we have sequenced the com- 
plete genomes of symbionts of the Aster Leafhopper (ALF), 
Macrosteles quadrilineatus (Cicadellidae: Deltocephalinae). 
This species is a widespread pest of North American agriculture 
(Frost et al. 201 1). As found in other members of Deltocepha- 
linae, it harbors Sulcia and the betaproteobacterial co-primary 
symbiont, Nasuia (Noda et al. 2012; Koga et al. 2013). 

Materials and Methods 

Material Acquisition and Sequencing 

Samples of M. quadrilineatus were field-collected from Yale 
West Campus, West Haven, CT, USA, on Symphyotrichum 
lanceolatum (Asteraceae). Specimens were identified accord- 
ing to Kwon (1988) and with the barcoding mitochondrial 
locus Cytochrome Oxidase I (following Le Roux and Rubinoff 
2009). In order to confirm the bacteriome association of the 
targeted bacterial lineages, whole-specimen Fluorescent 
in situ Hybridization was performed following Koga et al. 
(2009). Fluorescing probes were designed as follows: Bet940 
(5 / -TTAATCCACATCATCCACCG-3 / ) for the Betaproteobac- 
teria labeled with 5' AlexaFlour-555 probe modification, and 
Sulc664R (5 / -CCMCACATTCCAGMTACTCC-3 / ) for Sulcia, 
labeled with a 5' AlexaFlour-647 modification (Invitrogen). 
Briefly, specimens were dissected in 70% EtOH, fixed in Car- 
noy's solution, and incubated in alcoholic hydrogen peroxide 
for 2 weeks. Material was then rinsed and stored in 100% 
ethanol, rehydrated in PBSTx, and rinsed with a hybridization 
buffer. Specimens were incubated in a hybridization-probe 
solution, containing the oligo-probes and DAPI counterstain. 
Specimens were imaged with a Nikon Eclipse TE2000-U epi- 
fluorescence microscope (fig. 2). 

Ten individual insects were pooled for genomic DNA 
(gDNA) sequencing. Total gDNA was extracted using Qiagen 
DNeasy extraction kit (Qiagen Corp.) in a 50jil volume. 
Genomic sequencing was done using lllumina's MiSeq plat- 
form. gDNA was sheared into library insert sizes of 500 bp, 
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Fig. 1. — Schematic summary of the evolution of symbiotic associations in major lineages of Auchenorrhyncha, with emphasis on the Cicadomorpha 
(cicadas, spittlebugs, leafhoppers, and treehoppers). The focal symbionts, Nasuia and Sulcia from Macrosteles quadrilineatus (Cicadellidae: Deltocephalinae), 
are demarcated with a blue and red arrow, respectively. Phylogenies show evolutionary events in symbioses based on molecular phylogenetic studies (Moran 
et al. 2003, 2005; Takiya et al. 2006; Urban and Cryan 2012; Noda et al. 2012; Bressan and Mulligan 2013; Hou et al. 2013; Koga et al. 2013). Events are 
color-coded as follows: Host insect phylogeny (from Cryan and Urban 2012) gray background; Sulcia, red; betaproteobacterial symbiont lineage 
{BetaSymb = Zinderia + Nasuia + Vidania), blue; symbiont loss, black dashed line; and symbiont replacement, green. Identified symbiont names are given 
at the tips. This schematic is not a full account of host relationships and symbiont associations, as many species and genera remain to be explored, and some 
inferred losses may reflect incomplete knowledge. Sulcia is hypothesized to have been associated with the Auchenorrhyncha since its emergence 260- 
280 Ma (see time scale and corresponding blue box; Moran et al. 2005). Although further evidence is required to test the hypothesis that the 
Auchenorrhyncha share a common betaproteobacterial symbiont (BetaSymb) that was acquired early in its diversification, this potential relationship is 
shown here with a blue dashed line. The hypothesis that the Cicadomorpha do share the closely related BetaSymb (Zinderia + Nasuia), as evidenced by our 
genomic and phylogenetic results, is shown in a solid blue line (as is Vidania found throughout the Fulgoroidea based on Urban and Cryan 2012). 



and library amplification (eight polymerase chain reaction 
[PCR] cycles) was performed with high-fidelity polymerase 
(KAPA Bio HiFi) to reduce known problems of biased amplifi- 
cation in AT-rich symbiont genomes (Aird et al. 201 1; Sloan 
and Moran 2012a). The prepared library was sequenced using 
paired-end orientation on a 2 x 250 bp lane. Library prepara- 
tion and sequencing were done at the Yale Center for 
Genome Analysis. 

Assembly 

Because sequenced samples include a complex metagenomic 
mixture of host and bacterial symbionts, we implemented a 
multi-tiered approach to produce preliminary genome contigs, 
collect associated bacterial reads, and reassemble and improve 



scaffolds (Sloan et al. 201 3). A first pass assembly was done on 
1/1 6th of the read data (1,937,157 reads) using Velvet 
vl.2.08 with a kmer length of 63 (Zerbino and Birney 
2008). Whole genome assembly returned nearly 12,000 con- 
tigs greater than 500 bp, the majority of which had relatively 
low kmer coverage. Velvet yielded 1 0 large contigs of over 
1 0 kb in length, which were identified using BlastN searches 
against the NCBI nr database. The three largest contigs, one at 
1 90 kb and two at -56 kb, were unambiguously assigned to 
Sulcia and to Betaproteobacteria, respectively (supplementary 
fig S1, Supplementary Material online). We hereafter refer to 
these symbiont lineages as Sulcia-AU 1 (of the Aster 
Leafhopper) and Nasuia-AU 1 . The remaining large contigs rep- 
resent insect nuclear DNA and the complete mitochondrial 
genome. 
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Fig. 2. — Fluorescence in situ hybridization of the bacteriome of 
Macrosteles quadrilineatus (Cicadellidae: Deltocephalinae), housing 
Sulcia muelleri (green; Bacteroidetes), and Nasuia deltocephalinicola (red; 
Betaproteobacteria). Blue fluorescence shows host nuclear DNA. (a) 
Localization of the bacteriome along the anterolateral abdominal seg- 
ments, (b) The dissected bacteriome comprising two distinct bacteriocyte 
types that house Su/c/a-ALF and Nasuia-AU 1 . White scale bars equal 
100>m (a) and 50 (im. (b) Note. — Red autofluorescence of the host tho- 
racic plates. 

To remove symbiont-derived reads from the total read 
pool, data were mapped to the Velvet contigs using SOAP2 
V2.21 (settings: -m 10 -x1000 -g 3 -r 2; Li et al. 2009). 
Isolated symbiont reads were then reassembled using MIRA 
(Chevreux et al. 1999) with paired-end insert size set to 
100-600 bp and normal quality settings. MIRA confirmed 
the 1 90 kb scaffold of Sulcia-AU 1 , but broke Nasuia-AU se- 
quences into three large contigs: 12, 42, and -54 kb. MIRA 
also produced a 5 kb contig including a 1 .1 kb tandem repeat 
that assembled with the ends of the two largest contigs, over- 
lapping each end by over 450 bp. The Velvet and Mira contigs 
were assembled into a consensus scaffold that filled most 
ambiguous sections introduced by Velvet (e.g., strings of 
Ns). This merged the two largest Nasuia-AU 1 contigs, leaving 
two large scaffolds of -102 kb and -12 kb in length. These 
were assembled using PCR and Sanger sequencing (primers: 
98836F-ACGCCGCCGAGAATAAGTG and 100182R-ATACAC 
AACGCGCCCTTCG). The remaining ends of the Nasuia-AU 
scaffold contained a 400 bp overlap that was further sup- 
ported by paired-end read data that bridged the two ends. 

To assess the quality of the merged scaffolds, to improve 
genome assemblies, and to determine overall coverage for the 
symbiont genomes, total read data were mapped to the sym- 
biont scaffolds using SOAP2. Read mapping was run itera- 
tively until it converged on a consensus genome map, using 
custom Perl scripts (Sloan and Moran 2012a). This corrected 
indels and base pair substitutions within initial assemblies. To 
verify complete coverage of the circular bacterial genomes, 
scaffolds were broken, reoriented, and reads were remapped 
to confirm closure. Sulcia-AU read coverage was relatively 



uniform, whereas Nasuia-AU mapping revealed a coverage 
spike of roughly 3-fold, starting at the nucleotide position 
55,964, and supporting the potential repeat region assembled 
by Mira. The 1.1 kb repeat is too large to verify with 
the 250 bp reads, so PCR amplification and Sanger se- 
quencing of the region was used to confirm the boundaries 
where the assembly initially broke (primers: 55201 F-ATTCGCC 
TCATCCCCTCTCG + 56233R-AGCCATAACCTTTTGCGTCGT) 
and to confirm the existence of the internal repeat with in- 
verse primers (primers: 56954F-AAATTAAGAGAGCTAGG 
+ 56087R-TTTAACATCTATCGCTTG). The copy number of ge- 
nomic repeats is difficult to verify and likely varies within insect 
hosts and within the symbionts that are polyploid (see Woyke 
et al. 201 0), as has been shown in Portiera, the endosymbiont 
of whiteflies (Sloan and Moran 2013). Thus, the reported 
genome includes a single copy of the repeating motif flagged 
with a repeat annotation. 

To determine whether other symbionts are present in the 
samples, and to search for potentially unintegrated Sulcia-AU 
and Nasuia-AU sequences, Velvet contigs longer than 500 bp 
were identified using translated Blastx searches against the 
NCBI nr database and taxonomically binned with MEGAN4 
(Huson et al. 2011). A vast majority of the reads were 
placed in insect groups or were unassignable; these likely rep- 
resent high copy insect genomic regions. Twenty-four contigs 
(500-700 bp in length) were close matches to regions of the 
genome of Arsenophonus nasoniae (Gammaproteobacteria). 
Arsenophonus is a cluster of bacterial symbionts present in 
many insect groups and in some plants (Novakova et al. 
2009). Given the extremely fragmented genomic contigs for 
Arsenophonus and its low relative coverage, it is unlikely to 
represent an obligate primary symbiont of Macrosteles. 

Genome Annotation 

First pass annotations were done with RAST and the Joint 
Genome Institute IMG ER pipelines (Aziz et al. 2008; Marko- 
witz et al. 2009), using the general bacterial translation table. 
These pipelines produced a relatively complete annotated 
Sulcia-AU genome sequence, but yielded many unidentifiable 
short open reading frames (ORFs) for Nasuia-AU. Previous 
research indicated that the small genomes of Hodgkinia and 
Zinderia use an alternative genetic code that has reassigned 
the UGA stop codon to tryptophan (McCutcheon et al. 
2009b; McCutcheon and Moran 2010). To see if this is the 
case for Nasuia, Glimmer3 was used to predict coding regions 
with the modified translation table 4 (Delcher et al. 2007). Use 
of the alternative table extended many previously truncated 
protein-coding genes to their expected lengths, providing 
strong evidence that Nasuia uses a UGA Stop -> Trp recoding. 
Many genes utilize UGA-Trp codons at UGG-Trp conserved in 
other species with the standard bacterial code (e.g., hisl, rpoB, 
dnaN, etc.; see also McCutcheon and Moran 2010). All in- 
ferred genes for both symbionts were manually checked and 
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curated using phmmer and hmmerscan searches with avail- 
able databases (NCBI nr, TIGRfam, SwissProt, and Pfam) in 
HMMER3 (Finn et al. 201 1). The EcoCyc and KEGG databases 
were used further to determine gene pathways and function 
(Moriya et al. 2007; Keseler et al. 201 1). RNA coding regions 
were identified using RNAammer, tRNAscan-SE, and searches 
against the Rfam database (Schattner et al. 2005; Lagesen 
et al. 2007; Burge et al. 2012). Genomes were manually an- 
notated with Geneious Pro v5.5.8 (Drummond et al. 2010). 

Phylogenomics 

To infer phylogenetic relationships of Nasuia-AU, we used 
Phyla_Amphora2, which identifies and aligns conserved geno- 
mic markers from a set of phylum-specific genes (Wang and 
Wu 2013). Conserved orthologs were determined using 
Amphora's marker scan with HMMER3 (E-value cutoff of 
0.01), which generated a 42 gene amino acid alignment 
(approx. 10k columns) for 102 Betaproteobacterial taxa (sup- 
plementary table S1, Supplementary Material online). The 
Betaproteobacterial symbiont of the mealybug Planococcus 
citri, Tremblaya princeps (placed in Betaproteobacteria: 
Burkholderiaceae; see Alves et al. 2013), was not included 
in the analysis as many of its genes were missing and our 
match threshold excluded nearly half of the available genes. 
Four outgroups were selected from Gammaproteobacteria 
and Alphaproteobacteria, and homologous genes were ex- 
tracted using Blastp searches. Individual gene sets were 
aligned using MAFFT v7.017 with the L-INS-i model (Katoh 
and Standley 2013), and manually curated to remove large 
gaps and ambiguously aligned regions. Amino acid likelihood 
substitution models were determined for each gene matrix 
with ProtTest3 (Darriba et al. 201 1). Phylogenetic trees were 
inferred using maximum likelihood criteria with RAxML v7.4.4 
(Stamatakis 2006). A concatenated matrix was partitioned by 
loci and assigned gene-specific GTR-based amino acid model. 
Tree searches were conducted for 1000 bootstraps with a 
final maximum likelihood search. 

Results 

Genome Assembly 

The lllumina MiSeq run yielded 30,994,516 quality-filtered 
reads. Sequence assembly and PC R verification produced com- 
plete genomes for Sulcia-AU and Nasuia-AU, with each 
genome comprising roughly 1.3% of total sequencing data. 
Coverage mapping of the merged assemblies for Sulcia-AU 
and the Nasuia-AU 1 produced unbroken scaffolds with high 
mean sequencing depths of 500 x and 850 x, respectively. 
Genome sizes for Sulcia-AU 1 and Nasuia-AU are 190,733 
and 112,091 bp (counting only a single copy of the 1.1 kb 
tandem repeat; see Materials and Methods), respectively. 
Both are the smallest in their phyla, and Nasuia-AU represents 
the smallest bacterial genome reported to date (fig. 3; 



McCutcheon and Moran 201 2). The GC content of both sym- 
bionts is low: Sulcia-AU at 24.0% and Nasuia-AU at 1 7.1 %. 



Sulcia: Gene Loss in a Stable Symbiont 

Sulcia-AU contains a genome of only 1 90 kb, whereas 
previous Sulcia genomes, from several xylem-feeding auche- 
norrhynchan species, are 245-280 kb in size (fig. 4; McCutch- 
eon and Moran 2007, 2010; McCutcheon et al. 2009a). The 
genome is perfectly syntenic with those of other Sulcia line- 
ages; however, protein-coding content is dramatically reduced 
(figs. 4 and 5). Sulcia-AU encodes 190 predicted protein- 
coding genes (five hypothetical and two pseudogenes), a 
single ribosomal RNA operon, noncoding RNA (rnpB), ssra 
RNA, and 30 tRNAs capable of recognizing all codons (fig. 5 
and supplementary table S2, Supplementary Material online). 
In contrast, Sulcia of the sharpshooter leaf hoppers and cicadas 
contain 226-228 and 246 genes, respectively (figs. 4 and 5; 
McCutcheon and Moran 2007; McCutcheon et al. 2009a; 
Woyke et al. 2010). Compared with previously sequenced 
Sulcia genomes, the Sulcia-AU genome retains the same set 
of orthologs for the synthesis of eight essential amino acids: 
leucine, isoleucine, threonine, lysine, arginine, tryptophan, 
phenylalanine, and valine (fig. 5) (McCutcheon and Moran 
2007, 2010; McCutcheon et al. 2009a; Woyke et al. 2010). 

Sulcia-AU has lost many genes underlying oxidative phos- 
phorylation (OXPHOS) and DNA replication and repair (figs. 4 
and 5). It lacks nearly all OXPHOS-related genes, including 
those encoding the complete cytochrome C oxidase complex, 
NADH dehydrogenase, menaquinone and ubiquinone biosyn- 
thesis proteins, and the molybdopterin oxioreductase complex 
that is involved in the terminal electron acceptor, DMSO re- 
ductase (Kisker et al. 1999). Sulcia-AU does, however, retain 
genes encoding F-Type ATP synthase machinery that could 
function to synthesize ATP from the built-up electrochemical 
gradient generated by missing enzymes. In the absence of 
most protein complexes in the OXPHOS pathway, its function 
remains unclear. 

Further gene losses involve DNA replication and repair 
(figs. 4 and 5), including genes encoding the DNA polymerase 
III beta (DnaN) subunit and DNA gyrase (GyrAB), which are 
responsible for the efficient replication of DNA by clamping 
the polymerase and relaxing coiling ahead of the DNA helicase 
(DnaB), respectively (Johanson and McHenry 1980; Sugino 
et al. 1980). Based on studies in model systems, these losses 
are expected to reduce replication efficiency (Burgers and 
Kornberg 1982; Smelkova and Marians 2001). Additionally, 
Sulcia-AU retains only a single gene (mutS) encoding a sub- 
unit of the methyl-directed DNA mismatch repair complex, 
which is responsible for correcting mismatches and frame- 
shift associated insertions and deletions in Escherchia coli 
(Schofield and Hsieh 2003; Polosina et al. 2009). The gene 
mutL is pseudogenized and is identifiable in Sulcia-AU as a 
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Fig. 3. — Genome map of the co-primary symbionts of Macrosteles quadrilineatus: Sulcia-AU (left) and Nasuia-AU 1 (right). Genomes are scaled relative to 
nucleotide count. The outermost track (blue) shows genes coded on the plus strand, whereas the adjacent tracks (teal) show genes coded on the minus 
strand. The other tracks show RNA genes, including the rRNA operons (purple), ncRNA (red), and tRNA (orange). The two graphs show GC skew (outermost: 
green = values above genome average, purple = below genome average) and GC content (innermost: gray = values above genome average, black = below 
genome average). Plots were created with DNAPIotter (Carver et al. 2009). 
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Fig. 4. — Synteny plot showing gene order conservation and gene loss between Sulcia-AU 1 (Macrosteles quadrilineatus) and 5c//c/a-GWSS (Homalodisca 
coagulata). Axes are in the number of nucleotides. Lettered breaks correspond to gene maps on the right hand side. Plots a-c show specific gene loss 
between 5ulcia-G\NSS (top) and Sulcia-AU 1 (bottom). Conserved orthologs are colored blue connected with gray shading, whereas genes that are lost in 
Sulcia- ALF are shown in red. 
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Fig. 5. — Genomic content of the co-primary symbionts from the leafhoppers, Macrosteles quadrilineatus (Sulda-AVF and Nasuia-AU 1 ) and Homalodisca 
coagulata (Sa/c/a-GWSS and Baumannia); and, the spittlebug, Clastoptera arizonana (Su/c/a-CARI & Zinderia). The major pathways include DNA replication, 
DNA repair, translation and processing, protein stability, and energy metabolism. Colored boxes indicate genes that are present, and white boxes indicate 
gene absence. Striped boxes indicate genes that are present, but possibly pseudogenized, or with uncertain function. 



135 bp segment, whereas homologs in other Sulcia lineages 
are 1 .1 kb in length (McCutcheon et al. 2009a). 

Nasuia: The Functions of the Smallest Symbiont Genome 

The Nasuia-AU 1 genome is 1 12 kb in size, 27 kb smaller than 
the next smallest known bacterial genome, that of Tremblaya 
princeps, a symbiont of the mealybug P. citri (139 kb; fig. 3; 
McCutcheon and von Dohlen 201 1), and 30 kb smaller than 
the genome of Hodgkinia cicadicola, the symbiont that co- 
resides with Sulcia in cicadas (143 kb; McCutcheon et al. 



2009b). The Nasuia-ALF genome encodes 137 predicted pro- 
tein-coding genes, a single ribosomal operon, and 29 tRNAs 
predicted to be capable of recognizing all necessary codons 
(fig. 5 and supplementary table S3, Supplementary Material 
online). Fifteen predicted ORFs are unassignable to class or 
function. Nasuia-ALF contains a 1.1 kb tandem repeat that 
encodes part of MetF, a hypothetical protein, and a cold 
shock protein (CspA). Similar large genomic repeats (ranging 
from 150 bp to 6.1 kb) and major structural inversions have 
been found in other reduced symbiont genomes including 
those of Carsonella (psyllids), Portiera (whiteflies), and 
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Tremblaya (McCutcheon and von Dohlen 2011; Sloan and 
Moran 2012b, 2013). 

Nasuia-ALF uses an alternative genetic code in which UGA 
is reassigned from Stop to Trp. Roughly 51% of all protein- 
coding genes contain at least one UGA codon. This reassign- 
ment also has been found in two other bacterial symbionts, 
Zinderia and Hodgkinia (McCutcheon et al. 2009b; 
McCutcheon and Moran 2010). Nasuia-AIF is further missing 
prfB encoding the peptide release RF2 that is responsible for 
recognizing the UGA stop codon (Capecchi 1967), as also 
observed in Zinderia and Hodgkinia. 

Nasuia-AIF retains the pathways for synthesizing both 
histidine and methionine, as found in other symbionts that 
co-reside with Sulcia (McCutcheon and Moran 2010). 
Nasuia-AlFs histidine pathway is complete. For methionine 
biosynthesis, Nasuia-AIF uses a cobalamin-independent 
direct sulfhydrylation pathway as previously found for 
Zinderia (McCutcheon and Moran 2010). Both retain CysHI 
enzymes for the sulfhydrylation pathway and methionine 
synthase (MetE; Hacham et al. 2003; McCutcheon and 
Moran 2010). Although both symbionts retain genes encod- 
ing homoserine acetyltransferase (MetX), they appear to use 
two different sulfhydrylases: Nasuia-AIF encodes succinylho- 
moserine sulfhydrylase ortholog (MetB), whereas Zinderia uses 
acetyl-homoserine ortholog (MetY; McCutcheon and Moran 
2010). Nasuia-AIF appears able to initiate acetyl-driven me- 
thionine biosynthesis initiated by MetX, but not the succinyl 
alternative. MetB is a versatile enzyme shown to utilize multi- 
ple homoserine-esterified intermediates, including both succi- 
nylhomoserine and acetyl homoserine (Hacham et al. 2003), 
and it may function to catabolize the latter. The plant-sap diet 
seems unlikely to furnish the missing metabolites, as plants use 
a phosphorylhomoserine pathway (Macnicol et al. 1981), and 
it has been suggested that the requisite homoserine is 
obtained from Sulcia's threonine synthesis pathways 
(McCutcheon and Moran 2007). 

Nasuia-Al? retains few other large pathways, but does 
retain genes underlying basic processes of replication, tran- 
scription, and translation, including genes encoding DNA 
and RNA polymerase subunits, RNA modification enzymes, 
ribosomal proteins, and some tRNA synthetases (fig. 5). It ap- 
pears able to perform a limited part of the OXPHOS pathway. 
It encodes the machinery to synthesize and assemble a cyto- 
chrome C oxidase, including both precursors (e.g., CcsB) and 
assembly proteins (e.g., SC01 -like proteins), although many 
of these proteins have incurred AT-biased substitutions and 
are difficult to identify specifically, or may be pseudogenes. 
Finally, Nasuia-AIF retains genes encoding several enzymes 
that function in the metabolism of carbohydrate and sugars, 
including flavodoxin reductase (Fpr) and phosphoglycerate 
mutase (GpmA); however, these pathways are greatly re- 
duced and appear incomplete. Other gene losses have left 
Nasuia without a complete set of aminoacyl tRNA synthetases 
{Nasuia-AIF has lost 14 additional genes for tRNA synthetases 



relative to Zinderia), most TCA cycle-related genes (e.g., sucAB 
and Ipd), and genes for several components of the DNA 
replication machinery {dnaG, dnaX, and holB). 

Phylogenomics 

Phylogenetic analyses based on amino acid sequences of 42 
proteins place Zinderia and Nasuia as sister lineages with high 
bootstrap support (BS = 1 00; fig. 6). Both lineages exhibit high 
levels of sequence divergence relative to other species of Beta- 
proteobacteria. Both are placed within the family Oxalobac- 
teriaceae (BS = 75). This placement was supported in previous 
phylogenetic results for Betaproteobacteria, including analyses 
based on 16S rRNA including Zinderia and Nasuia representa- 
tives (Noda et al. 2012; Koga et al. 2013) and analyses based 
on proteins from genome sequences including only Zinderia 
(McCutcheon and Moran 2010; Alves et al. 2013). 

Discussion 

The Origins of the Auchenorrhyncha Symbionts 

Macrosteles quadrilineatus is associated with two obligate 
symbionts, 5i//c/a-ALF and Nasuia-AU 1 , as previously observed 
for some other members of the Deltocephalinae (Noda et al. 
201 2; Ishii et al. 201 3). Most groups of Auchenorrhyncha are 
associated with Sulcia; however, symbionts that co-reside with 
Sulcia vary, implying multiple evolutionary transitions in major 
auchenorrhynchan clades (see fig. 1; Moran et al. 2003; 
McCutcheon et al. 2009b; Koga et al. 2013). Phylogenomic 
results and other biological characteristics support a close re- 
lationship of Nasuia with Zinderia, which co-resides with Sulcia 
in most cercopoid spittlebugs (Koga et al. 2013). Sequence- 
based phylogenetic analyses should be interpreted cautiously 
in the inference of bacterial symbiont relationships, as appar- 
ent relationships might reflect from a combination of conver- 
gence and long branch attraction due to the extreme AT 
nucleotide bias, rapid sequence evolution, and ancient diver- 
gence of symbionts relative to other bacteria (see McCutch- 
eon and Moran 2012). Despite this, our genomic results 
provide additional biological support for a close relationship 
between Nasuia and Zinderia. First, Nasuia and Zinderia share 
cobalamin-independent, direct methionine synthesis path- 
ways not found in the other co-primary symbionts of the 
Cicadomorpha (McCutcheon and Moran 2010). Second, 
and more distinctively, the reassignment of UGA observed in 
Nasuia and Zinderia has occurred only a few times in the his- 
tory of life, including the small AT-rich genomes of mycoplas- 
mas, several mitochondrial lineages, and Hodgkinia (Yamao 
et al. 1985; Knight et al. 2001; McCutcheon et al. 2009b; 
McCutcheon and Moran 2010). The reassignment is absent 
from most other reduced, often AT-rich genomes, such as 
those of Blattabacterium, Blochmannia, Buchnera, Carsonella, 
Portiera, Uzinura, and Tremblaya (Shigenobu et al. 2000; Clark 
et al. 2001 ; Gil et al. 2003; McCutcheon et al. 2009b; Sabree 
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Fig. 6. — Relationships between Nasuia deltocephalinicola and other Betaproteobacteria. Phylogenetic tree was inferred with maximum likelihood criteria 
for 107 taxa for 42 proteins (approx. 10,000 amino acid columns). Circles at nodes indicate bootstrap support values: black = >75, gray =55-75, and 
white = <55. See supplementary table S1, Supplementary Material online, for complete taxon sampling, and supplementary figure S2, Supplementary 
Material online, for full phylogeny. 



et al. 2009, 2013; McCutcheon and von Dohlen 201 1; Sloan 
and Moran 2012b). 

Several hypotheses exist to explain the evolution of codon 
reassignment, including the complete extinction of the UGA 
stop codon and its reemergence as a Trp codon, along with 
UGG (Osawa and Jukes 1988). However, even small and 
AT-biased genomes contain many UGA Stop codons. An al- 
ternative hypothesis, proposed for Hodgkinia, is that the reas- 
signment is triggered by an inactivating mutation in RF2 (prfB), 
which causes some peptides to be extended, but is not lethal 
(McCutcheon et al. 2009b). Such a mutation would resemble 
other deleterious mutations and gene inactivations that are 
fixed through genetic drift in reduced symbiont genomes. 
Selection would then favor replacing UGA stop codons with 
fully functional stop codons (UAA or UAG), freeing up UGA to 
bind with the tRNA-Trp anticodon. The tRNA-Trp may evolve 
to enhance the weak binding already present between its 
anticodon and UGA. Because all cases of this recoding occur 



in genomes with extensive gene loss (Nasuia, Zinderia, 
Hodgkinia, Mycoplasma, some mitochondrial genomes), the 
initial loss of prfB as a triggering event is a potential 
commonality. 

Our results bolster support for a symbiont clade containing 
Nasuia and Zinderia, and suggest that this BetaSymb clade 
represents an ancient associate of Sulcia, with these two sym- 
bionts co-residing in hosts at least since the origin of the 
Cicadomorpha (figs. 1 and 6). Preliminary phylogenetic evi- 
dence further suggests that the fulgorid symbiont, Vidania, 
also may be derived from a common ancestor shared with 
Nasuia and Zinderia (Urban and Cryan 2012; Koga et al. 
2013). Both Vidania and Zinderia have co-diversified with 
their respective host groups for at least 200 million years 
(McCutcheon and Moran 2010; Urban and Cryan 2012). 
We hypothesize that Vidania and Zinderia represent subclades 
of BetaSymb, and share an ancestor that was present in the 
shared ancestor of the fulgoroid and cicadomorph lineages 
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prior to their split 260 Ma (see fig. 1; Scherbakov 2002). A 
genomic sequence of Vidania, not yet available, could provide 
further evidence for or against this hypothesis. 

In earlier work on Auchenorrhyncha symbionts, based lar- 
gely on light microscopy, Buchner (1965) and his student 
Muller (1940, 1962) hypothesized the existence of the symbi- 
ont later characterized with molecular data and named Sulcia 
(Moran et al. 2005). They further hypothesized companion 
symbionts for different auchenorrhynchan groups, but, lack- 
ing molecular evidence, did not recognize common ancestry 
of betaproteobacterial symbionts present in several auchenor- 
rhynchan groups. 

The hypothesis that Sulcia and BetaSymb were co-resident 
in the shared ancestor of Cicadomorpha implies that replace- 
ments of BetaSymb with alternative symbionts have occurred 
multiple times and at different evolutionary depths (summa- 
rized in fig. 1). These transitions include switches to Hodgkinia 
(Alphaproteobacteria) in cicadas and to the Sodalis-Wke sym- 
biont (Gammaproteobacteria) in philaenine spittlebugs (Koga 
et al. 2013). Within the leafhoppers, the replacement of 
Nasuia with Baumannia appears to be correlated with the 
origin and diversification of the monophyletic xylem-feeding 
sharpshooters, 25-40 Ma (Moran et al. 2003; Takiya et al. 
2006). Baumannia retains one of the largest obligate symbiont 
genomes to date, encoding nearly complete pathways for 
energy metabolism and vitamin biosynthesis (Wu et al. 
2006), and contrasting with the reduced genome of Nasuia- 
ALF in other leafhoppers. Xylem feeding creates serious chal- 
lenges, requiring physical and metabolic compensations to 
overcome negative tension and low energy content 
(Brodbeck et al. 1993; Zimmermann et al. 1994; Novotny 
and Wilson 1997). By providing new gene repertoires, the 
acquisition of a novel bacterial symbiont may have enabled 
hosts to overcome these challenges, allowing the diversifica- 
tion of xylem-specialized insects (Takiya et al. 2006). However, 
symbiont replacements and losses have sometimes also oc- 
curred in phloem-feeding lineages. For example, Scaphoideus 
titanus, another phloem feeder classified in the Deltocephali- 
nae, lacks both Sulcia and Nasuia and instead possesses trans- 
ovarially transmitted Cardinium (Bacteroidetes) symbionts plus 
a yeast-like symbiont related to those found in certain lineages 
of planthoppers and aphids (Sacchi et al. 2008). Some other 
leafhopper groups, such as the subfamily Typhlocybinae that 
feeds on more nutritious parenchyma plant tissues, appear to 
have lost their symbionts entirely (Buchner 1 965; Moran et al. 
2005). Nutritional aspects of these reported cases of loss or 
replacements have not been investigated. 

Extreme Conservation of Essential Amino Acid Synthesis 

A main role of these bacterial symbionts is to provision host 
insects with the essential amino acids that their hosts can 
neither synthesize, nor obtain in sufficient quantities from 
their plant sap diets. In general, plant sap offers low 



concentrations of the essential amino acids relative to the 
nonessential amino acids (Brodbeck et al. 1993; Sandstrom 
and Moran 1999). Essential amino acid synthesis by the sym- 
bionts of M. quadrilineatus is achieved by the conserved com- 
plementation pattern previously observed in dual symbiont 
systems: Sulcia-AIF synthesizes 8 of the 10 required by 
the host, whereas Nasuia-Al? synthesizes the remaining 2 
(reviewed by McCutcheon and Moran 2010). The only excep- 
tion to this pattern reported so far is in spittlebugs, where 
Sulcia-CAR\ has additionally lost the tryptophan pathway, 
which is retained by Zinderia along with histidine and methi- 
onine synthesis pathways (McCutcheon and Moran 2010). All 
previous examples of this complementation in the Auchenor- 
rhyncha were from xylem feeders (McCutcheon and Moran 

2010) . Our study shows that a similar pattern occurs in 
phloem feeders. Although phloem is far more nutritious 
than is xylem, having 10- to 100-fold higher concentration 
of free amino acids, the amino acid profile is usually unbal- 
anced relative to animal nutritional needs, and it varies with 
plant species and growth stage (Sandstrom and Moran 1 999). 
A high relative concentration of single amino acids can render 
sap noxious or unusable to insects (Brodbeck and Strong 
1987), and insects have been shown to perform better on 
sap in which amino acids' profiles resemble experimentally 
determined nutritional requirements for insects (Dadd 1985). 

The amino acid pathways in the symbionts of M. quadrili- 
neatus, and also in the other Auchenorrhyncha, are relatively 
complete and do not appear to require the host or companion 
symbionts to contribute intermediate steps. The only pro- 
posed exception is the transfer of homoserine produced by 
Sulcia to its co-primary symbionts for biosynthesis of methio- 
nine (McCutcheon and Moran 2012). In contrast, in symbiont 
systems in the sister suborder, Sternorrhyncha (e.g., psyllids, 
whiteflies, aphids, and scale insects), some pathways are miss- 
ing intermediate and terminal catabolic steps. For example, 
many primary symbionts in the suborder are missing the 
aspartate amino transferase (AspC) required for the final 
step in phenylalanine synthesis, and, with the exception of 
Carsonella in psyllids, they have also lost the branched-chain 
aminotransferase (HvE) required to complete the synthesis of 
isoleucine, leucine, and valine (Shigenobu et al. 2000; 
Nakabachi et al. 2006; McCutcheon and von Dohlen 2011; 
Sloan and Moran 2012b). For the few systems for which 
genome sequences are available for insect hosts, host en- 
zymes capable of completing these pathways are encoded 
(Wilson et al. 2010; McCutcheon and von Dohlen 2011), 
and these genes are upregulated in aphid bacteriocytes 
(Hansen and Moran 201 1; Poliakov et al. 201 1). In addition, 
the synthesis of tryptophan and/or arginine is accomplished by 
pathway integration between co-occurring symbionts in psyl- 
lids (Sloan and Moran 2013), the aphid Cinara cedri (Gosalbes 
et al. 2008), and mealybugs (McCutcheon and von Dohlen 

2011) . The fact that all sequenced symbionts of the 
Auchenorrhyncha, despite millions of years of divergence 
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and several replacements, encode relatively autonomous 
pathways suggests that the host and companion symbionts 
cannot contribute interspersed steps in the synthesis of essen- 
tial amino acids. The complementary, but separate, synthesis 
of amino acids in the Auchenorrhyncha may reflect the an- 
cient segregation of symbionts into different bacteriocytes 
(Buchner 1965; Moran et al. 2005; Koga et al. 2013). 

Finally, the difference in the methionine synthesis pathways 
between Zinderia and Nasuia remains difficult to interpret, as 
only these two genomes are available for betaproteobacterial 
symbionts of Auchenorrhyncha and amino acid sequence 
identity for both loci to available sequences is low. 
(Zinderia's MetY is 55% identical to that of Herbaspirillum 
and Nasuia's MetB is 35% identical to that of Burkholderia.) 
Given the hypothesized antiquity of Zinderia and Nasuia, it is 
possible that their shared ancestor had both copies and that 
one or the other has been lost in the different lineages. Several 
available genomes for Betaproteobacteria indicate that multi- 
ple copies of sulfhydrylase can be present in members of 
Oxalobacteraceae (family in which Zinderia and Nasuia are 
placed) (Pedrosa et al. 201 1) and the related Burkholderiaceae 
(Holden et al. 2004). An alternative hypothesis, that one of the 
sulfhydrase genes was acquired through horizontal transfer, 
seems less likely, as reduced symbiont genomes generally do 
not acquire foreign genes and the gene sequences display no 
signs of foreign origin. 

The Shared Loss of Energy Metabolism Genes 

The most striking loss of genes from Sulda-AIF and Nasuia- 
ALF involve the pathways underlying energy metabolism and 
particularly the OXPHOS pathway. Both symbionts retain only 
a single OXPHOS gene complex, and neither retain genes re- 
quired to synthesize NADH dehydrogenase, which is a unique 
loss in M. quadrilineatus' symbionts compared with symbionts 
of previously studied auchenorrhynchan species. These losses 
indicate that both symbionts have a reduced capacity to syn- 
thesize ATP. Although ATP demands may be reduced, as basic 
cellular functions have been stripped from both symbionts, 
ATP is required for numerous pathways in both symbionts, 
raising the question of how they generate energy. 

Speculatively, the loss of these functions, relative to other 
auchenorrhynchan symbiont systems, is related to the en- 
riched sugars and energy-rich nature of phloem. Carbohy- 
drate concentrations in phloem are several orders of 
magnitude higher (~1000x) than those in xylem (Mattson 
1980). Evidence from the excreta of the xylem- and phloem- 
feeding insects indicates that less than 1 % of xylem feeders' 
waste is organic carbon, but that it may be as high as 70% in 
the waste of phloem feeders (van Hook et al. 1 980; Andersen 
et al. 1989). This excess of sugar in phloem sap potentially 
reduces the need for OXPHOS-generated ATP in symbionts, 
leading to the gross reduction of these pathways in the 
Sulda-AIF and Nasuia-Al? genomes. As expected under this 



hypothesis, the genomes of the other xylem-feeding Auche- 
norrhyncha do retain more complete OXPHOS pathways, in- 
cluding the missing NADH dehydrogenase (Wu et al. 2006; 
McCutcheon and Moran 2007, 2010). 

In the Sternorrhyncha, OXPHOS genes are variably retained 
in symbiont genomes. Buchnera of aphids and Portiera of 
whiteflies have near-complete sets of OXPHOS-related 
genes, but several have lost the ATP synthase complex (e.g., 
Shigenobu et al. 2000; Perez-Brocal et al. 2006). Similar losses 
are observed in the Tremblaya-Moranella dual symbionts of 
the mealybug P. citri (McCutcheon and von Dohlen 2011). 
Sequenced genomes of Carsonella from several psyllid species 
reveal that they have completely lost NADH and other genes 
from OXPHOS protein complexes (Sloan and Moran 2012b). 
Sets of genes lost show some overlap between the symbiont 
systems of phloem-feeding Auchenorrhyncha and Sternor- 
rhyncha, although differences are evident. Perhaps much of 
energy metabolism is a nonessential symbiont trait that the 
host is able to supplant when pathways are lost. 

In order to metabolize phloem's carbohydrate-rich diet, 
symbionts would need to retain the genetic mechanisms to 
operate the glycolytic, pyruvate metabolic, and citric acid 
pathways. Both Sulda-ALF and Nasuia-AIF do possess some 
genes for metabolizing carbohydrates; and, for Sulcia-CAR\ of 
cicadas, several of these (e.g., gapA, Ipd, sucAB) are among 
the most highly expressed genes based on proteome analyses 
(McCutcheon et al. 2009a). However, these genetic reper- 
toires in Macrosteles' symbionts are incomplete, calling into 
question whether or how well each symbiont can metabolize 
sugars. Possibly, the host exchanges intermediate metabolites 
required for carbohydrate metabolism directly with the bac- 
teriocytes. Insects can metabolize sugars independently of 
their symbionts, and they maintain carbohydrate and fatty 
acid stores in nearby fat tissues (Kunieda et al. 2006; Arrese 
and Soulages 2010). In the obligate intracellular symbionts of 
the grain weevil, Sitophilus zeamis (Curculionidae), host- 
encoded sugar transport proteins are indeed expressed in 
bacteriocyte tissues (Heddi et al. 2005). It is hypothesized 
that these transporters allow sugars to enter bacteriocytes 
along a chemiosmotic gradient where they are metabolized 
in response to the carbohydrate-rich nature of the host in- 
sect's diet (e.g., wheat; Heddi et al. 2005). 

An alternative explanation is that the host may supply its 
symbionts with ATP directly, or that the required missing 
OXPHOS genes are provided by one of the symbiotic partners. 
Suida-AlV retains F-type ATP synthase, which requires an elec- 
trochemical gradient to drive the synthesis of ATP (Kasimoglu 
et al. 1 996). However, in the absence of all other members of 
the OXPHOS chain, it is not clear how the electrochemical 
gradient is achieved by Suida-AlV. Possibly, the existing 
OXPHOS genes are shared between co-primary symbionts 
because Suida-AlV and Nasuia-AlV encode complementary 
components, although how these would be integrated is 
not clear. Proteins or protein complexes may be shared 
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between the host and symbiont partners through the endo- 
plasmic reticulum (ER), which is closely associated with host- 
derived membranes that encase intracellular symbionts (re- 
viewed by Houk and Griffiths 1980). The loss and retention 
of single, yet complementary, genes in the OXPHOS pathway 
by Sulcia-AIF and Nasuia-AIF is a mystery, requiring further 
work to understand symbiont and host contributions to the 
metabolism of energy in this symbiosis. 

Supplementary Material 

Supplementary figure S1 and tables S1-S3 are available at 
Genome Biology and Evolution online (http://www.gbe. 
oxfordjournals.org/). 
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